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Abstract 

We calculate the average multiplicity of jets in e^e~ annihilation as a 
function of both the jet resolution parameter, j/cut, and the thrust T. Our 
result resums to all orders the leading and next-to-leading logarithms in ?/cut 
and 1 — T, and is exact up to second order in as. This allows a comparison 
between the number of jets found using jet algorithms and the ability to 
distinguish different jet topologies via thrust. 
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1 Introduction 



The study of jets and their physics in e"'"e~ annihilation has added much to our under- 
standing of perturbative QCD. Although hadrons are the final state particles seen in 
such collisions, it is quarks and gluons whose dynamics are described by perturbative 
QCD. The definition of jets, via infrared and collinear safe jet clustering algorithms, 
bridges the gap between the theoretically accessible partons and the experimentally 
observed hadrons. 

The multiplicity of jets in e"'"e~ annihilation events provides a good example of this. 
By contrasting the jet multiplicity with the hadron multiplicity we may see some of the 
advantages of using jets as our object of study. In particular, when predicting the number 
of hadrons in the final state, one must invoke some non-perturbative hadronization model 
at an unknown hadronization scale Qq. This results in an arbitrary overall normalisation, 
so that only the energy variation of the hadron multiplicity can be calculated within 
perturbative QCD. In contrast, by examining the jet multiplicity, the hadronization 
scale is replaced by the known jet resolution scale Qo = Q^/y^, and no arbitrary 
parameters are introduced. Therefore, the jet multiplicity is fully calculable within 
perturbative QCD, including its absolute value. In fact, since the jet resolution scale 
is firmly within our control it can be varied at will, effectively studying the energy 
dependence of the multiplicity in a single experiment. 

For large values of the jet resolution scale, the jet multiplicity is reliably predicted by 
fixed-order perturbation theory. However, the appearance of logarithms of ?/cut spoils the 
perturbative expansion for finer jet resolutions and these logarithms must be summed to 
all orders. With the advent of the kx_ (or Durham) jet finding algorithm it is possible 
to sum the leading and next-to-leading logarithms of ?/cut to all orders, making the jet 
multiplicity also reliable for small resolution scalesp|. This is the jet algorithm that we 
will use throughout this paper, and it defines jets according to the following iterative 
procedure: 

• For each pair of particles calculate a separation defined by. 



where Ei and Ej are the particles' energies and Oij is the angle between their 



• For the pair of particles with the smallest separation, yij, combine their momenta. 
Pi and pj, together to form a pseudo-particle of momentum pij. In the ii^-scheme, 
Pij = Pi+Pj- 

• Repeat the above two steps until all the separations of particles or pseudo-particles 
are less than the jet resolution parameter ?/cut- The remaining particles and pseudo- 
particles are then called jets. 

This algorithm can be applied to form jets from either partons (for theoretical consider- 
ations) or hadrons (for experimental analysis). The difference between performing the 
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momenta. 
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algorithm at the partonic and hadronic levels is expected to be suppressed by at least 
one power oil/Q . 

The theoretical calculation and experimental analysis of event shapes have also con- 
tributed greatly to our understanding of perturbative QCD. In particular, the thrust 
and heavy jet mass of an event have been extensively studied[^ and show good agree- 
ment with experiment 0. It is interesting to calculate the multiplicity of jets in e^e~ 
annihilation, retaining the jet kinematics so that the multiplicity can be expressed for 
specific values of one of these event shapes. In this paper we calculate the multiplicity 
of jets in e"'"e~ annihilation as a function of the thrust, summing all large logarithms to 
next-to-leading logarithmic accuracy. This was done for the multiplicity of hadrons in 
Ref. 0. The thrust for an event is defined by, 

T = max/5^^^|, (2) 

where n is a unit vector chosen to maximise the right-hand-side, and Pi are the three- 
momenta of the final state partons or hadrons. Of course, our results can be easily 
applied to other similar observables, such as the heavy jet mass. 

The structure of the paper will be as follows. In section 2 we shall consider the 
fixed-order contribution to the multiplicity up to C(ag). Then, in section 3, we will 
resum large logarithms of r = 1 — T and ?/cut to all orders in as, and to next-to- leading 
logarithmic accuracy, using the coherent branching formalism. To avoid double counting, 
we will discuss the matching of this result to fixed order in section 4, and finally, we 
present numerical results in section 5. 



2 Fixed Order 

We begin by defining the quantity N{r, i/cut) as the multiplicity of subjets resolved in 
the k± algorithm at a scale j/cut, in events with thrust greater than 1 — r. It is generally 
given by, 

A^(r,ycut) = f """^^'^^"'^ (3) 

n=2 ^tot 

where cr„(r, ycut) is the cross section for events with thrust greater than 1 — r and exactly 
n jets, and atot is the total cross section for the production of hadrons. We will require 
this up to second order in the strong coupling constant, as. Defining (T^™')(r, ?/cut) to be 
the 0{a^) contribution to crn(r, ?/cut), the multiplicity to 0{al) is given by, 

N{T,ycnt) = 

2 (q-Q + (r, ycut) + (^2^ (r, ^cut)) + 3 {a^^^ (r, ^cut) + 0-3^^ i^, Vent)) + 4 ^ (r, ^cut) 

(4) 
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where (Jq is the Born cross-section and C is some number that we will not eventually 
require. Truncating at 0{a'^), we have: 

. as ^A"sV , fi 1 (1). X , 1 (2). ^ , fas\ 
1 C [ — +1 —a, ' r,u„itJ H cto ' r,w„,t, + — 



N{t, ycnt) = 1 --O'l — I + li -I — a2 ' Kj, ycut j H 0-2 ' l^, l/cut J . , , 

TT \7r/\7r/cro ctq Vtt/ 

+ 3 ((1 - ^) -ai'HT,|/cut) + -ot\T,y.uS\ + 4-af(T, I/cut). (5) 

V V TT / (7o (To / (To 

Since at C((3;s) every event has either 2 or 3 jets and at 0{(x^ every event has 2, 3 or 
4 jets, we can write, 

(7W(r) = a«(r,|/,„0+4'^(r,l/cut), (6) 
= (r,|/cut)+af (r,|/eut)+af (r,|/cut). (7) 

This allows us to write the multiplicity in terms of more inclusive quantities. Further 
simplification is made by rewriting our expression in terms of the cross section for events 
with thrust less than 1 — r, cT(r), using, 

a«(r) + a«(r)=(^)(To, (8) 

(T(^)(r) + a(^)(r)=(^)'cao. (9) 

This gives, 

A^(r,ycut) = 2(l-(l-^)laW(r)-la(^)(r)) 

V V TT / (Jq (Tq / 

VV TT / (To (To J Gq 

= 2(l-i?(T))+i?3(T,?/cut)+2i?4(T,?/cut), (10) 

where R{r) is the ratio of the cross section for events with thrust less than 1 — r to the 
total cross section, and Rnir, ycut) is the ratio for events with n jets and thrust greater 
than 1 — T, evaluated to the appropriate order in q;s, i.e. 

-R(t) = — O-(t), -Rn(T,ycut) = — 0-n(T,t/cut)- (H) 
Ctot CTtot 



^.i Zeroth Order 

At zeroth order, 0(Q;g), i.e. with no gluon emission, the cross section is given by the 

production of a quark-antiquark pair, emitted back-to-back in the centre of mass frame. 
The quark and antiquark will always be placed in separate jets, so the jet multiplicity 
will be two, and the thrust axis will fall along their common axis, trivially giving a 
thrust of one. Therefore we have, 

7V(r,|/cut) = 2. (12) 
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2.2 First Order 

We easily obtain the multiplicity to 0{as) by truncating Eq. (pUj), giving, 

(13) 

Since neither of these cross sections involve the singular region of thrust equal to 1, they 
can be calculated in four dimensions. At this order, i.e. 0{as), the only contribution is 
from the partonic process e"'"e~ — > qqg. 

For -R(t), the cross section for producing events with thrust greater than 1 — r divided 
by the total cross section to hadrons, one obtains, 

R(r) = ^{21ogV + 31ogr + 5_!^ 

- 6r log r - 4 log(l - r) log r - 6r - 3(1 - 2r) log(l - 2r) 

+ 4Li2 (^) - + 21og^(l - r)|, (14) 

where only the first line contributes as r ^ 0. 

Next we turn to i?3(r, ?/cut), the cross section for producing three jets with thrust 
greater than 1 — r divided by the total cross section to hadrons. To this order, 

R3{r,ycut) = — 4^V>l/cut). (15) 

The contributing partonic process is still e"'"e~ — >■ qqg, but now we demand that the 
parton separations be greater than ?/cut in order to maintain a three jet final state. 

For ycnt < 1=7 we have, 

i?3(^,Z/cut) = — — < / dxi / dX2 

Zir [Ji-T Ji-xi/2 



'l-2/join j (1 - Xi)(l - X2)(Xi + X2 - 1) 

= RDiVcut) - Rir), (16) 

where, 

\/ycut(8 + ?/cut) - Vcut ,^ ^, 

2/join = ^ , (17) 

and -R_D(z/cut) is the three-jet fraction in the algorithm at a scale ?/cut, given in Ref. 0. 
For ?/cut > this becomes, 

„ , . Cpas f^~y^-^ , f ^\7ylu'~^f , 2(x? + xi + (2 - xi - xs)^^) 
R3{r,ycut) = — — / dxi / dx2- 



2tT Ji-t Jl-xi/2 (1 — Xi)(l — X2)(Xi + X2 — 1) 



271 Jl-T 1 — Xi 
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with 

gixuVcnt) = -{2- 3xi(l - X,)) log ^^""^^y"^;;^^ - 3(2 - x.f | " ^'-"^ ~ (19) 

Xl Vent - [l-Xiy 1 + ycnt - 

This integral has a closed analytic form, but does not give a very compact expression, 
so we do not reproduce it here. 

In order to match our resummed result to fixed order we must also extract the 
logarithmic behaviour for r, ?/cut 1- We obtain: 



Cpas r 2 log' ^ - 3 log + log 64 7/eut > 2t' 



2n \ log' ^ - 2 log' i - 3 log ^ + log 64 + log^ 2 y,,t < 2r^. 

(20) 

Note that these two solutions are continuous at the matching point, and that it can be 
shifted, with logarithmic accuracy, to: 

f 21og'^-31og^ + log64 > r^ 

- ^ I log' ^ - 2 log' i - 3 log ^ + log 64 yeut < r', ^^1) 

or: 

Rsir, ycut) = ^ (2 log' ^ - 3 log — - e(r' - ^eut) log' — | . (22) 

2n [ Vent Vent Vent J 

2.3 Second Order 

Since we could make the first-order calculation in four dimensions, it is possible to 
calculate the second-order terms using a standard Monte Carlo NLO program like 
EVENT2||]. 

Now Eq. (^) requires no modification and the multiplicity is, 

iV(r, y,,t) = 2 (1 - Rir)) + R^{t, y,^,) + 2R,{t, y,^,), (23) 

where -R(r) and -R3(r, ?/cut) should now be evaluated to NLO, and -R4(r, ?/cut) to leading 
order. 

Notice that at this order we can study the dependence on the renormalization scale, 
since we have 



with b = — |A^/, and likewise for Rsi^r, y cut)- 
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3 All-Orders Summation of Large Logarithms 

In this section we present the all-orders resummation of all leading (ag log™" r log'""™" ?/cut) 
and next-to-leading (ag log"" r log'"""""^ 2/cut and log™"V log'''"™ |/cut) logarithms, 
which appear for small ycut and/or r. These logarithms can be summed to all orders 
using the coherent branching formalism. 
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To next-to-leading logarithmic accuracy, we have, 



(25) 



Here, -Fg(/c^, Q"^) is the probability that a quark formed at scale Q has a mass below k, 
and log Fq has been calculated in Ref . to next-to-leading logarithmic accuracy. Here 
we only require the expression for F„ itself to next-to-leading logarithmic accuracy, which 



IS, 



where. 



«s(g), 



2A) log 



2(1 -A) log 



1 - 2A 



1 - A 



/2(A) 



A 



(26) 

(27) 
(28) 



Also, Nq{k'^,Q'^;QQ) is the contribution to the multiplicity in a quark jet formed at a 
scale Q and resolved at a scale Qq, from jet masses below k. 

In order to derive Ng^k"^, Q^; Qq), let us first consider the quantity n^^k^, Qq) dk^ 
which is the multiplicity of partons of species a in the jet, where the squared mass of 
the jet lies between k'^ and /c^ + dk^ . Clearly, 



(29) 



a={q,9} ' 

This obeys an exclusive evolution equation given by. 



<(fc', Q'; Ql) = a(q2) 5,. 6{e) e(Q - Qo) 

+eiQ-Qo)l -^^Idz PUz)J^ dkgj^ dk, 

{n-^ikl z\^)f,{kl (1 - zfq') e{zq - Qo) 

+nl{kl (1 - zfq'- Ql)f,{kl z\')Q{{l - z)q - Qo)) 

5(e-z{l-z)q'-^-^\. (30) 

This has a simple physical interpretation. Firstly, it is possible that the quark does 
not emit any radiation. Then only the original quark will contribute to the multiplicity 
if it is sufficiently hard (i.e. Q > Qo). Clearly the resulting jet will be massless. This 
gives rise to the first term. Alternatively, the quark may survive until it reaches a 
scale q before emitting a gluon of mass kg and light-cone momentum fraction z (and 
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leaving the quark with mass kg and hght-cone momentum fraction 1 — z). Then the 
resolvable partons originating from both quark and gluon must be counted. The function 
fq{k'^, Q^) dk"^ in the above is the probability that a quark formed at a scale Q has a 
squared mass between k'^ and k"^ + dk'^, i.e. it is the derivative of Fq{k'^, Q'^) with respect 
to k"^. A is the Sudakov form factor, 



logA(Q=) = - /; ^ ( dzP„{z r'^':_''' . (31) 



v2 ^1 tts(^(l - z)q) 

Jo g2 

These equations have only formal meaning, because the integrals are divergent and 
the form factor is zero. To make them physically meaningful, they can be calculated 
explicitly by imposing an infrared cutoff z{l — z)q > e, or by working in c? = 4 — 2e 
dimensions. They are however infrared finite, so the limit e ^ can be taken smoothly. 
Thus we always imply such a procedure. 

Differentiating with respect to logQ^ yields an inclusive evolution equation: 



Ql) ^ j; ..^^^Mi^P,.(.) {/; dki f 

(n^(fc2, z'Q'; QDUkl (1 - zfQ^)Q{z^Q^ - Ql) 



(1 - zfQ'- Ql)f,{kl z^Q')Q{{l - zfQ' - Q\ 



S - Z{1 - Z)Q' - f - ^ j - Q2. g2)@(g _ g^)| _ (33) 

The mass of the gluon will not significantly contribute to the mass of the jet so to next- 
to-leading logarithmic accuracy it can be omitted from the 5-function. This 5-function 
and the normalisation conditions, 
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defq{k\Q^) = 1, (33) 
dk\{k\Q';Ql) = X^{Q';Ql), (34) 

can then be used to integrate out fc^ and k"^. Here J\fg{Q'^;Ql) is the multiplicity of 
partons of species a inside a quark jet, where the quark is emitted at a scale Q and the 
partons are resolved at a scale Qo and no demand is made on the mass of the jet. This 
has been derived in Ref. 0, 

Afg^{Q';Ql) = M-iQ';Ql) + l^^mQ';Ql), (35) 

^f^{Q';Ql) = ^fHQ';Ql)-l^^m';Ql), (36) 



where. 



A^+(g'; Qo) = ^1 (j^y [Ib+M)Kb{zo) + Kb+i{z,)Ib{zo)] 



Zo\;iCA b 



^f-{Q'■,Ql) 

AfiQ';Ql) = (-) [lB{z^)KBizo)-KB{z,)lB{zo)]. (37) 

.Zi , 



6 V 3 3 3Ca J 

ly and Ky are the modified Bessel functions and we have used the notation of Ref. [Q] 
where, 

2 _ ?,2txCa 2 _ 327rC^ 

Removal of these integrations gives, 



+n^((l - z){e - z{l - z)Q^), (1 - ^)2g2. g2)@((i _ ^)g _ g„)0(A;2 _ _ ^)g2) 

-n^(fc2,Q2.g2)0(g_g^))_ (40) 

The above evolution is dominated by the emission of soft gluons (z ~ 0), so to next-to- 
leading logarithmic accuracy z can be replaced by zero in any smooth functions. Since 
neither fg nor ra^ are smooth functions (they contain 1 / fc^-like terms) such a replacement 
cannot be made in their arguments. However, their integrated distributions, 

,2 ^2\ / T„2 f /2 ^2\ 



F,ie,Q') = / dq'U{q\Q'), (41) 

•J 

N^g{k\Q'-Ql) = dq'n-^{q',Q';Ql), (42) 

which are the functions we are really interested in, are smooth and obey the same 
evolution equations as fg and n° but with different boundary conditions. 

Integrating Eq. (^Oj) over k"^ and replacing z with zero where appropriate gives to 
the required accuracy. 



U^{z^Q^; Ql)Fg{e - zQ^, Q^)Q{zQ - Qo)e{k^ - zQ^) 

+N^{e - zQ^ g2)0(^2 _ ^g2)@(g _ g^) 
-N^{k\Q';Ql)e{Q-Qo)). (43) 

Furthermore, k"^ — zQ^ can be replaced with fc^ in Fg and A^^ to next-to-leading loga- 
rithmic accuracy, giving the final evolution equation in integral form. 



^1 asizq) 
uz 

q" Jo 



+ r'^r^^^^..W^.(^^^')-^;(^V;Qg)e(fc^- V)e(.g-go) 



fc2 q^ Jk'yq^ 2n 



2 ,^,dz^^Pggiz)N^ik\Q';Qf,). (44) 
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This has a formal solution: 
N^{k',Q';Ql)=expi^- 



as{zq) 
fc2 Jk^/q^ 2tc 



+ 



dq^ 

— exp 





dz^^^P,,{z)\[N^^{e-,Ql) 

F,{k\q') 



'''4f dz^-^P,,iz 

fc2 q^ Jk^/q2 



vr 



X dz'^^Pg,{z)M^{z\^-Ql)Q{e - zq^)Q{zq-Q,)]. (45) 
The integrals are considerably simplified by inserting the unintegrated expression for Fq [^], 



Fq{k\Q')=e^v[- r % t dz 

[ Jfc2 q^ Jk^/q^ 
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(46) 



(47) 



so that the expression for A*"" becomes, 

N^q{k\ Q': Ql) = Fq{k\ Q') { M^q{k^-, Ql) + Cl{k\ Q'- Ql)} , 
where, 

q{k\Q';Ql)= r% ['dz^^P,,iz)Af^izY;Ql)Qie-zq')Qizq-Q,). (48) 
Jfc2 q jQ ztt 

The remaining integral can be done by making a change of variable from to 
q^ = z'^q^. To next-to-leading logarithmic accuracy, 

Cl{k\Q^-Ql) = t'"^' '4 dz^^P,,{z)M^{q'-,Ql)Q{q-Q,) 

Ji) q J q/Q ^71" 



= 2Ci. 
+ 2Cf f 

Jk 



q^ 27r 
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'k' 



Af^ie;Ql)iogi-Je{q-Qo). (49) 



Making the B-functions more explicit and rearranging gives. 



■Kir-^Ql) log 



Q 



+ Q{e -QI)2Cf 
- e(P - QQo)2C^ / 



2 g2 27r ^ 



Ae(g';Qg)log 



g2 27r 

C^{r(A;2;Q2)_ja(^4/g2.g2)|^ 



A/';(g';Q^)iogU 



(50) 
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where 



Ql) = Q{e - Ql) fl Ql) log 



(51) 



This integral is most easily done using the change of variables, Eq. (^9[), with the 
obvious addition, 

2 327rC^ 



so that it becomes. 

Using the results, 

'''"^^U^{z-z,){zl-z^) - 

20 Z 

dz 



1 



dz 



2Ca J Z 



mz;zo){zi-z') 



2(^+{zk;zo)-l)+ABAr{zk;zo 



Z 



Af{z; Zo){zi - z') = 2M{zk;zQ)~zi/zl + l 



(52) 
(53) 

(54) 
(55) 



we obtain, 



0(^fc - 2:0)7^ \M+{zk] zq)-1 + 2BJ\f{zk] zq] 



C 



2Af{zk]Zo)- zllzl + l 
C 



nk\Ql) = e{z,-zo)^{2Af{zk;zo)-zl/zl + l) 



where the parameters B and C are given by Eq. ([ 

3 b Ca 



and 



(56) 
(57) 

(58) 



However, the quantity of interest is the multiplicity of all jets, independently of their 
flavour, i.e. = + N^. Therefore we finally have, 

N,{k\ Q'; Ql) = F,{k\ Q') { X,ik'; Ql) + {l{k\ Ql) - lik'/Q', Ql)) } , (59) 



with, 



Hk',Ql) 



e{Zk - Zq)^ [Af+iZk] ^0) - 1 + '^BAfiZk, Zq) 

+ (5-1) {2Af{zk;zo)-zl/zl + l] 



(60) 



The above equation conforms with naive expectations. One might expect that the 
number of jets of mass below k"^ resulting from a quark created at a scale Q^ would be 
simply the probability of finding a jet of mass below fc^ multiplied by the number of jets 
within. Indeed, this is the first term of our expression for A^^^, and the naive expectation 
requires only the addition of next-to-leading logarithmic corrections. 
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For completeness, we also present the result for Ng{k^,Q^,Ql), the multiphcity of 
partons found within a gluon of squared mass below k"^. The derivation closely fol- 
lows that of the quark case above and will not be reproduced here. For the tagged 
multiplicities we have, 

iV;(A;^ Q'- Ql) = F,{k\ Q') Ql) + Ca {r{k\ Ql) - r{k^/Q\ Ql))} , (61) 

Here Fg{k'^, Q"^) is the integrated gluon jet mass distribution [Q, and takes the same form 
as Fq, Eq. but with the coefficients of /i and /2 now being Ca and 6/2 respectively. 
Also AfgikP'] Qq) is the multiphcity of parton species a found within the gluon0. 

Notice that the function is identical to that in the quark case. This is because the 
only contribution to these functions is from the singular parts of the splitting kernels in 
the evolution equation. The appropriate kernel for the quark is Pgq{z) = 2Cf/z + 0{1) 
and for the gluon is Pgg{z) = 2Ca/z + C(l) (the other kernels Pgg and Pqg have no 
singular parts and cannot contribute). The only difference between the two cases is the 
colour factor, Ca or Cp, which is reflected in the coefficient of 

Clearly, the untagged multiplicity is given by, 

Ngik^ Q'; Ql) = Fg{k\ Q') { Mg{k'-, Ql) + Ca {l{k\ Ql) - I{k'/Q\ Ql)) } . (62) 
4 Matching to Fixed Order 

Of course, the resummed result also contains part of the fixed-order contribution, which 
has already been included. It must be matched to fixed order so that we do not over- 
count. This is done by expanding the resummed result to next-to- leading order in as, 
and subtracting the offending piece from our result. Expanding each term individually, 
we obtain: 

i^.(e-V,Q^) ^ 1 - g (C^L^ - f CpL) 

+ (g)' [\ClL' - [\hCF + lCl]L' + ilhCF + lCl]L^) , (63) 

where as = (^s{Q)i 

J\fq{e-e-'k') ^ i + ^(^ICfI'~ICfI) + (^^J (^l^CFCAl'-^CF[3CA-b]P 
2{Ca - Cp)CFTnNf{8CFTnNf - 2CATnNf - Cj) 

where as = as{k), and finally, 

iie;e-^e)^^{\i^) 

^/«s\Vi^ ,4 , i,;3 2(C^ - CF)T^Nf{llCl - ACpTnNf) 
where again as = as(^)- 
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Note that, in the above, the coefficients of al-^^ and a'^P are not beheved to be 
correct, since they are only next-to- next-to-leading logarithms. Nevertheless, they must 
be considered when performing the matching since the fixed-order result is exact to 
order ag. It only remains to write the various as values in terms of as{Q)'- 



«s(fcvg) 

27r 



as 
2n 



2%) 



1 1 

) logT^ 
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2tx) 



log 



fc2 



(66) 
(67) 



The result of this expansion is most easily written in the form of Eqs. ( [T0| ) and (|2 
with the additional notation: 



asiQ) 
2% 



+ 



2tx 



(68) 



Then, 



A{e- 



. ^-l-L\ 



9(^2 T 



2CpL^ — 3CpL, 

-2ClL^ + Cf{QCf + h)L^ - Cf{^ + \Cf)L^ 
{2CfI^ - 3CFl)eil) - CFil - LfQ{l - L), 

\hL - (2CfL^ - SCfL) + ICf){2CfI^ - 3CfI) 

+^CfCaI'^ + i^b — ^Ca)CfI^ 
8{Ca - CF)CFTRNf{6Cl + CaTrNj - QCFTRNf)f 

{hL - {2CfL^ - 3CfL) + ICf)Cf{1 - if 



(69) 
(70) 
(71) 



0(0 



+^^CfCa{1 - Lf + \hCF{l - Lf 
A{Ca - CF)CFTRNf{llCl - 4CFTnNf){l - Lf 

9Cl 



Q{l-L). 



(72) 



These have been checked by solving the evolution equation iteratively, verifying that the 
claimed solution does actually satisfy the evolution equation, at least to second order 
in as- 



5 Numerical Results 

Combining the fixed-order and resummed results, we obtain the contribution to the 
mean jet multiplicity from events with thrust between T and T + dT, n(l — T, ycut)dT. 
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1-T 

Figure 1: The mean number of jets in e^e~ annihilation as a function of thrust, using 
the 0{a\) result alone (dashed) and matched with the resummed result (solid). 



To obtain the mean jet multiplicity as a function of thrust we normalise this to the 
number of events in the same range, r(l — T)dT, to give 

{N){T) = (73) 

shown in Fig. |1] for a representative value of |/cut = 10~^. We fix /i = y/s = niz and 
Q^s(^z) = 0.120. We see that for 1 — T ^ Hcut, the resummation is extremely important. 
For 1 — T ~ ?/cut, the physical threshold, the NLO corrections are negative and even 
after resumming large logarithms to all orders, give an unphysical result, with the mean 
number of jets falling below 2. This region is anyway outside perturbative control, since 
it corresponds to A;^ = s(l — T) ~ Ay/s. 

In Fig. 1^ we show the result for fixed thrust T = 0.95 as a function of ?/cut- The 
resummation is again seen to be important for ycut ^ 1 — T. 

The thrust can be used to separate two-jet from three-jet events. In Fig. ^ we show 
the multiplicity of jets in each sample. These are defined by 

(A^)(T>t) = ^^^~^^p\ (74) 

and 

mr < t) = ^''^"fr^'ll'j^)-^-' . (75) 

respectively. They can be directly compared with the results of Ref. where the 
separation into two- and three-jet events was made using the same jet algorithm as 
the one in which the jets were counted — the k± algorithm. Note that for small r/cnt 
the 0{a'^) results appear to be unphysical with the multiplicity in two-jet events being 
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Figure 2: The mean number of jets in e^e annihilation as a function of |/cut for fixed 
thrust, using the 0{al) result alone (dashed) and matched with the resummed result 
(solid). 




Ycut 

Figure 3: The mean number of jets in e~^e~ annihilation as a function of i/cut for events 
with thrust above and below a fixed cut, using the 0{al) result alone (dashed) and 
matched with the resummed result (solid). 
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Ycut 

Figure 4: The mean number of jets in e~^e~ annihilation as a function of ?/cut for fixed 
thrust, using the resummed result matched with the 0{as) (dashed) and 0(a;g) (solid) 
results, for renormalization scale choices /i = a/s/2, ^/s and 2^/s. 



higher than for three-jet events. This is somewhat misleading as they are effectively 
calculated to different orders owing to the different denominators in Eqs. (|7^) and ([75|) , 

The all-orders results are physically behaved, because the first uncalculated term in 
either result is suppressed by at least two orders of as and two powers of log(?/cut)- 

One might naively expect that the number of jets in a three-jet event be larger than 
in a two-jet event by a factor {2Cf + Ca)/{2Cf) ~ 2, which is certainly not the case 
in Fig. ^ where for small T/cut it is around 1.3. This is comparable to the ratio found in 
the kx algorithm of around 1.4 and, as discussed in Ref. is because the form factor 
suppression for the gluon jet is always stronger than for the quark jets. 

Finally, in Fig. ^ we show the variation of the result at fixed thrust under changes 
of renormalization scale between \/s/2 and 2y/s. We see that the 0{a\) result matched 
with the all-orders resummation is considerably more stable than the 0(as) plus re- 
summed result. 
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